Flame fronts in Supernovae la and their pulsational stability 
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The structure of the deflagration burning front in type la supernovae is considered. The parame- 
ters of the flame are obtained: its normal velocity and thickness. The results are in good agreement 
with previous work of different authors. After that the question of pulsational instability of the 
flame subject to plane perturbations is considered. The flame can be unstable if hydrodynamics can 
be ignored, e.g. in solid-body propellants. However, with account of hydrodynamics we find that 
the flame in type la supernovae is pulsationally stable with realistic parameters of reactions and 
thermal conduction. 



I. INTRODUCTION 

Supernovae explosions are among the most spectacular 
events in the Universe: their energy release significantly 
mixes the interstellar medium and acts like a driving force 
in gas dynamics of galaxies and production of cosmic 
rays. The luminosity of an exploding star becomes com- 
parable with the luminosity of the progenitor galaxy, and 
allows one to make observations of processes in the most 
distant regions of the Universe. If a regularity of param- 
eters from explosion to explosion is found in any subtype 
of supernovae, it could make possible to realize quantita- 
tive measurements of distances and to apply supernovae 
to cosmology. 

Despite the long history of investigations of these 
events the complete understanding of underlying phys- 
ical processes is still missing. There are several types of 
supernova explosions that very much differ physically and 
are related to different objects. Here we will consider only 
one subtype - thermonuclear explosions, the so-called su- 
pernovae of type la, SNIa. Such explosions are usually as- 
sociated with degenerate stars, white dwarfs. The study 
of SNIa shows that the underlying physical phenomenon 
creating the explosion is the thermonuclear burning of 
premixed carbon-oxygen fuel. 

The mode of the explosive nuclear burning in super- 
novae is still a controversial issue, in spite of many years 
of the research in the field. Four decades ago, Arnett 
1] was the first to model the supersonic combustion, 
i.e. detonation, in supernovae. Later, Ivanova et al. [2[ 
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obtained a sub-sonic flame (deflagration) propagating in 
a spontaneous regime with pulsations and a subsequent 
transition to detonation, and Nomoto et al. |3[ modelled 
the deflagration propagating due to a convective heat 
transfer. Both detonation and deflagration have their 
merits and problems in explaining the supernova phe- 
nomenon [see, e.g. 0j- It is not clear if detonation suc- 
ceeds or fails to develop, but it is clear that in any case the 
combustion must be much faster than suggested by the 
analysis of the propagation of a laminar one-dimensional 
flame. 

From the microscopic point of view the one- 
dimensional nuclear flame is a wave described essentially 
in the same way as it was done by Zeldovich and Frank- 
Kamenetsky [5| in spite of complications introduced by 
nuclear kinetics and a very high conductivity of dense 
presupernova matter. It is found that the conductive 
flame propagates in a presupernova with the speed which 
is too slow to explain the supernova outburst correctly: 
the flame Mach number is of the order of percent or less 

a. 

The fuel consumption can be naturally accelerated by 
the development of the instabilities inherent to the flame 
front. As it is explained in the classical paper by Landau 
0, HI , the hydrodynamic instability leads to wrinkling or 
roughening of the front surface, and hence to an increase 
of its area with respect to the smooth front and conse- 
quently to acceleration of the flame propagation. In ex- 
treme cases, the instabilities can lead to a transition from 
the regime of slow flame propagation to the regime of 
detonation. Since the flame propagates in gravitational 
field, and the burned ashes have lower density than the 
unburned fuel, the Rayleigh-Taylor (RT) instability is of- 
ten considered to be the dominant instability governing 
the corrugation of the front [i-El. The RT instability 
gives birth to a turbulent cascade providing for an accel- 
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eration of the flame front. However it leads to additional 
difficulties in modelingthe SN event [T3 - [l8| . 

It is well known 0, 0, [H, H3 that large portions of a 
slow planar flame front are unstable with respect to the 
large scale bending. This universal instability is called 
the Landau-Darrieus (LD) instability. For the wave- 
lengths much longer than the flame thickness it does 
not depend, on complex processes which take place in 
the burning zone. Development of the LD instability de- 
pends only on the sign of Ap = p u — pt>, where p u and 
Pb are densities of the unburned and burned "gases" re- 
spectively. The LD instability of planar flame fronts with 
respect to large scale bending takes place if, and only if 
Ap>0. 

The LD instability plays an important role in many 
physical phenomena such as the usual chemical burn- 
ing of gases, explosive boiling of liquids [2l|, electroweak 
phase transitions [22j , dynamics of thermally bistable ga s 
[23j , and thermonuclear burning in supernovae [2J, |25| . 
The detailed consideration of non-linear stage of LD in- 
stability and the calculation of the fractal dimension of 
the flame front for that case is given by Blinnikov and 
Sasorov [2(|, Joulin [27|. It is interesting, and somewhat 
puzzling, that a similar dependence of the flame fractal 
dimension on the density discontinuity was found in the 
3D SPH simulations of the flame subject to RT instability 

Both LD and RT instabilities develop on scales much 
larger than the flame thickness and they can be success- 
fully studied in the approximation of the discontinuous 
front. This approximation is not valid for another in- 
stability, first discovered by Zeldovich [29[ in his inves- 
tigation of the powder combustion. This instability is a 
result of a strong temperature dependence of the reac- 
tions rates, so that the local fluctuations of the heating 
rate caused by the temperature fluctuations cannot be 
controlled by the thermal conduction. It can lead to a 
pulsating regime of the front propagation and to a renor- 
malization of the mean front velocity [30] . This instabil- 
ity can develop even for one-dimensional perturbations 
when the plane front preserves its shape. Let us denote 
it as TP (thermal-pulsational) instability. After paper 
by [3l| it was studied quantitatively in many papers [see 
references in [32[ . 

The paper is organized as follows: in the Section |TT] 
the physical conditions and a model of WD explosion are 
presented, in the Section IIIII the stationary propagating 
flame is considered analytically and numerically. Some 
conclusions on the dependence of the results on nuclear 
network are presented. In Sections II VI and IVl the stability 
of the flame subject to plane perturbations is considered. 



II. THE MODEL 

A very nice review on the subject is given by [33[ [see 
also H3, [lH ■ For self contained presentation we will sum- 
marize some essential points here. There exist several 



scenarios of SNIa explosions. The most popular are: 
single-degenerate scenario, double-degenerate scenario, 
and sub-Chandrasekhar mass explosion. In this paper, as 
we think, the most reasonable one, the single-degenerate 
scenario, is considered. According to it a binary stellar 
system, which is a progenitor of the supernova, consists 
of a white dwarf (WD) and a non degenerate star. Dur- 
ing accretion of matter on the WD it approaches the 
Chandrasekhar mass limit and at some moment becomes 
unstable. In the language of equation of state it happens 
because the adiabatic exponent approaches the critical 
value 7 « 4/3. During this process the temperature in 
the centre of the WD rises and nuclear burning of mat- 
ter (that is degenerate 12 C and 16 O according to evo- 
lutionary models) is ignited. The temperature of igni- 
tion depends on matter density and can be found in [36| . 
For p ~ 10 9 g/cm 3 it is about T ~ 10 s K. But this 
burning is very slow and does not propagate outwards 
from the centre until its energy release is compensated 
for by various losses (the most significant are the neu- 
trino losses). The dynamical stage sets in later, when 
T rises up to 10 9 K and the flame is born. This is the 
beginning of the supernova explosion. The flame starts 
to propagate from the centre of the star to its surface. 
The regime of the flame propagation is under intensive 
investigation and the answer is not yet found. But it is 
not fully unknown, the observations imply some limita- 
tions on it. There exist two types of stationary regimes: 
1) deflagration when velocity is small compared to the 
speed of sound and the flame is driven by dissipative ef- 
fects: thermo-conductivity or diffusion; 2)detonation, i.e. 
a supersonic wave with the shock front where the tem- 
perature jumps up drastically leading to fast burning. If 
the burning of the whole star proceeds in the detonation 
regime, it burns up to Fe-peak elements, and that contra- 
dicts observations: in a real supernova about half of the 
star should consist of the intermediate elements. Pure de- 
flagration regime does not succeed too: the star expands 
with velocity, which is faster than the flame velocity so 
the temperature drastically drops down and all the burn- 
ing terminates. The only feasible successful regime is the 
mixed one: the flame starts with deflagration in high 
density matter, where the expansion coefficient is small. 
Then it somehow accelerates at the radius, which is usu- 
ally characterized by some critical density, and passes 
to detonation. This model theoretically succeeds in ex- 
plaining all the parameters of the explosion, but has one 
tuning parameter, critical density, p cr it- From compari- 
son of simulations with observations it is obtained that 
p cri t ~ 2 • 10 7 g/cm 3 [see, e.g. H, H3- But the global 
problem in SNIa physics (that is surely not solved in this 
paper) is to construct the model of the explosion with 
no tuning parameters. One should calculate the critical 
density of the burning regime transition from the first 
principles, this is the problem for the future. 

Let us discuss the physical conditions in the white 
dwarf. Suppose its mass is close to the Chandrasekhar 
mass limit. The density in the centre of WD is p ~ 10 9 
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g/cm 3 . The chemical composition is mainly 12 C and 
16 O. So the temperature at which the flame can burn 
under the specified conditions is about T ~ 10 9 K. Two 
main physical processes that matter in SNIa explosions 
are thermo-conductivity and nuclear reactions (burning) . 
In this case ions are non-degenerate and non-relativistic, 
electrons on the contrary are strongly or semi-degenerate 
and relativistic. EOS for such matter is presented in 38 1 
and takes into account ions, all degrees of degeneracy of 
electrons and photons. 

Under these conditions the thermo-conductivity plays 
the leading role in comparison with other dissipative pro- 
cesses. It consists of two components: radiative and 
electronic. The electronic thermo-conductivity, whch we 
use in this paper, is considered in (39l. |40|. An approx- 
imate equations for radiative thermo-conductivity are 
presented in 41[. We can estimate the value of the dif- 
fusive effects in the system. Despite the fact that in the 
white dwarf the matter is well mixed, let us consider the 
value of the diffusion coefficient. The ion diffusion can 
be approximated as 



D — XiV T - 



(1) 



Here vt is the thermal speed of ions, A,; is the ion mean 
free-path, which may be crudely estimated as the inter- 
ion distance. So the Lewis number is about Le = 10 4 . 
This result can be easily explained by the differences 
of speeds: relativistic particles contribute to thermo- 
conductivity, and non-relativistic - to diffusion. For vis- 
cosity 



r\ = mnXiVT- 



(2) 



So the Prandl number is Pr = 10 -4 . 

At the temperatures specified earlier the nuclear re- 
actions proceed in a branched network with a lot of iso- 
topes. For simplification let us suppose that WD consists 
only of 12 C. Then the first reaction in the nuclear net- 
work is 12 C+ 12 C— > 24 Mg* (here "*" means exited nuclei, 
it is unstable and decays in 3 channels with p, n or a 
in final state). Its caloricity is q = 5.6 • 10 17 erg/g, so 
it leads to significant temperature rise, up to T ~ 10 10 
K, and can provoke further burning. In this paper two 
variants of nuclear network are considered: 

1. A simplified network with only one reaction 
12 C+ 12 C^ 24 Mg* [the rate can be founding. Ac- 
cording to the fact that it is the first reaction in the 
network and electromagnetic nature of this reaction 
it can be a good approximation for the whole net- 
work at the first glance, and all dynamics will be 
determined by this reaction. The complete burning 
up to Ni could also be modeled in this framework 
by fixing the rate and changing the caloricity to 
q = 9.2 • 10 17 erg/g. 

2. a-chain with 13 isotopes from (43|, we will call it 

"APROX13". 



We will make also additional simulations with initial 
16 composition using APROX13. Because of the high 
matter density the degeneracy parameter satisfies T ~ 
-Seoul /i^kin ~ 1, so the electron nuclear screening effects 
should be taken into account. The accurate screening 
factor in a non-ideal gas can be obtained only by Monte 
Carlo simulations, and is presented, e.g., in [44]. We 
should emphasize also that many-order discrepancy in 
the screening factor can arise because of incorrect defini- 
tions, see details in [40| . 

Based on the above consideration we can write down 
the hydrodynamic system of equations which should be 
solved: 

dp 

— + div(pv) = 0, 



djfiCj) 
dt 

dpv 

~~dt 

dE 



- div(pvv) = -Vj> + /, (3) 
,-dbr(v(E + p) - kVT) = )>//,,/>',,. 

ij 

p = p{p,c l ,e). 

III. ONE-DIMENSIONAL FLAME PROPERTIES 
A. The Problem Set-up 

The main questions we pose here are as follows: What 
is the normal velocity of flame propagation, what is the 
expansion coefficient of matter over the flame front, and 
what are the physical parameters of matter after the 
flame traversal. To answer these questions we should 
study the plane burning front microphysics taking into 
account the processes on the scale of the thickness of the 
front. Let us make some theoretical estimates. For the 
stationary deflagration regime the characteristic times of 
burning and heat transfer over the front should be: 



Q _d 2 

Tnucl — • 3 7"cond — j 
S K 



(4) 



so the flame thickness and velocity can be estimated as: 




d 

T 




(5) 



In table U the estimated velocity values are presented 
for different conditions in WD, when the rate of burn- 
ing S is determined by one reaction proposed earlier, 
(2 12 C^ 24 Mg*). 

The exact front parameters can be obtained only by 
the numerical simulations of the system (|3]) in ID. For 
the general rates J2y of the reactions in system ([3]) the 
density and temperature distributions in a steady burn- 
ing wave are unknown, so the initial conditions with the 
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TABLE I: Theoretical estimate of the front velocity and thick- 
ness with 2 12 C-> 24 Mg* (Eqns. @-©) : 



p, g/cm 3 


v, km/s 


Ax, cm 


2 • 10 s 
7- 10 8 
2 • 10 9 


203 
742 
1300 


2.5 • 10" 5 
3.0 • 10" 6 
1.0 • 10~ 6 



steady flame could not be imposed. We will ignite the 
flame with a warm wall and try to eliminate all interfering 
perturbations. We expect that the stationary flame that 
will appear in such conditions is a natural phenomenon. 
The following task for numerical calculations is consid- 
ered: The region of interest is filled with a uniform mat- 
ter distribution. For simplicity it will be 12 C with density 
Po and temperature To, chosen so that the characteristic 
time of the flame burning with To is much larger than 
the dynamical time in the problem under consideration. 
In this case the ignition is controlled by the boundaries. 
The right hand side boundary condition sets a constant 
external pressure. The left hand side boundary condition 
is a hard wall v = 0, and its temperature rises with time 
according to the linear law: 

2"^ 2^0 

Teft wall = T H 1. (6) 

T 

The choice of r will be clarified below, Ti is the temper- 
ature which is larger than the temperature of the active 
burning. All these quantities depend on the task and are 
tuned "by hands" after inspection of the trial values in 
each case. Let L be the size of the region of interest (it is 
smaller than the whole region of calculation). The burn- 
ing leads to the temperature growth and therefore the 
pressure rises as well, which generates sound waves. For 
the sake of avoiding unnecessary perturbations we should 
give time for sound waves to quit the region of interest: 
r L/c s . In this case the pressure in the region L is 
constant with a good accuracy. 

Let us describe the numerical methods for the solutions 
of system ([3]). Physically, the system consists of 3 parts: 
hydrodynamics, thermo-conductivity, and burning. For 
the solution we will use the method of splitting in phys- 
ical processes. As an example consider the equation for 
energy E in system ([3]). At each time step the equation 
is split as: 





= div( K VT). (10) 

cond 



The equation for cj is split analogously (in other equa- 
tions of system © only one physical process in one 
equation is presented). Eqns. (|7l)-([T0l) are solved in 
3 sub-steps: The hydrodynamic part in which Ci — 
const and no energy generation is assumed; the nuclear 
part (p, v — const); and the thermo-conductivity part 
(p, Ci, v — const). 

For the solution of hyperbolic PDE's we will use two 
solvers (for cross-check): the implicit numerical scheme 
with Newton solver described in (4||, and Kurganov- 
Tadmor scheme [46[. Parabolic PDE (fT0|) is solved with 
Crank-Nicholson numerical scheme. The kinetic equa- 
tion for Ci with simplified nuclear network can be easily 
integrated for a small time step. The APROX13 nuclear 
network has its own integrator. 

Criteria that are used for determination of time and 
space discretization of parameters are the following: 

At = min(A£ courant ,i, Ai nucU ), (11) 

where Ai COU rant,i is the Courant condition for i-cell and 
Atnudj is the time step at which all reagent concentra- 
tions change not more than by 1%. 



B. Results of simulations 

Now we will present the results of our computations. 
The flame speed is determined by the time evolution of 
the coordinate of the centre of the front x{t) (by defi- 
nition it is the point, where cqi2 = 0.5). An example 
of the time dependence, x(t), is presented in Fig. [1] In 
Fig. [2] the coordinate dependence of the temperature for 
a sequence of times is shown. We see that some time af- 
ter the ignition the flame stabilizes and propagates with 
constant velocity. Fitting the x(t) dependence we obtain 
the velocity of flame relative to the wall, v. The normal 
flame velocity is then v n = ptvj p u . Table ILT1 presents the 
results of simulations for different parameters of initial 
matter conditions for all nuclear network variants. Com- 
paring our results with those of Q we can claim a good 
(but not excellent) agreement. 

Fig. [3J shows coordinate dependence of ci for different 
elements at fixed time for APROX13. The numerically 
obtained results for 1-step nuclear network are in good 
agreement with theoretical estimates (table fl} . The re- 
sults with APROX13 show that the flame speed is about 
an order of magnitude slower. It is an important result 
showing that the pure carbon burning to Mg is not a 
good approximation, despite the fact that it dominates 
at early stages of the burning (see Fig. [3]). 

Flame speed obtained with APROX13 (which is 
supposed to be closer to the real speed) can be 
calculated theoretically if we substitute the rate 
of reaction 24 Mg(a, 7) 28 Si in Eq. §5$ instead of 
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TABLE II: Results of the flame simulations: 



Compos. 


Po, 


g/cm 3 




-tmax, 1U 1^ 


Pu/Pb 


Vn, km/s 


Axfr, cm 




2 


ID 


i\ /r _ 

Mg 


0.9 


1 C /I 

1.54 


(U.l 


1. ( • 1U 








iNl 


/ .y 


1.00 


Izz 


l.z ■ 1U 








APROX13 


6.8 


1.60 


18.2 


3.0 • 10 -4 




r-7 
1 


ID 


i\ /r _ 

Mg 


n i 
9.1 


1.33 


302 


l.o • 1U 








iNl 




1.0 1 




1.0 ■ 1U 








A DD AV 1 Q 

ArKUAlu 


o.O 


1 ^ 

1.O0 


00.^ 


^± • _LU 




2 


10 9 


Mg 


11.3 


1.26 


854 


5.5 • 10 -6 








Ni 


13.8 


1.40 


1241 


1.0 ■ 10~ 5 








APROX13 


9.8 


1.23 


134 


1 ■ KT 4 


16Q 


2 


10 8 


APROX13 


6.1 


1.38 


0.94 


2 ■ 10 -2 




2 


10 9 


APROX13 


8.8 


1.17 


22.1 


5 ■ 10" 4 




0.0008 



FIG. 1: An example of the front coordinate dependence on FIG. 2: Sequential profiles of temperature for ti < ti < £3 < 
time for p = 2 • 10 9 g/cm 3 ; APROX13 nuclear network is used. t4, p = 1 ■ 10 9 g/cm 3 , APROX13 nuclear network 



2C— !>Mg. The speed should decrease by the factor 
^i?CM g (T 9 k, 10)/i? M gSi(79 « 10) ~ 10. This result is in 
accordance with simulations. Reactions with a particles 
dominate the whole burning and have more or less the 
same rates (the charge Z of heavy ions in the governing 
reactions does not vary significantly), so this estimation 
is accurate for all a-reactions. 



IV. THERMAL INSTABILITY OF NUCLEAR 
FLAMES 

Let now consider thermal-pulsational instability with 
switched-off hydrodynamics. In this case, because of ab- 
sence of sonic waves, this instability can be described 
purely without interference with other physical phenom- 
ena. Two conditions are necessary for the development 
of the TP instability [30I . l32j . First, thermal diffusion 
must be more than an order of magnitude larger than 
mass diffusion. This condition is undoubtedly satisfied 
in presupernovae where Le 3> 1. Second, the Zeldovich 



number [4j 



Ze = 



d In q 
<91nT 



(12) 



P, Xi 



which characterizes the temperature dependence of the 
heating rate q, must be high enough. Here P is the 
pressure, Xi (i — 1,2,...) are the abundances of the 
reactants, and the derivative is evaluated at the tem- 
perature of the burning matter. The thermal instability 
takes place when Ze is higher than a certain critical value 
Ze c . From the Arrhenius law, and in the approxima- 
tion when the reaction zone is assumed to be negligibly 
thin in comparison to the preheat zone of the flame, it 
is found Ze c = 4 + 2y/5 = 8.47 [13]. The recent numeri- 
cal value |48| , obtained by relaxing the approximation of 
delta function kinetics, is Ze c ~ 8.24 — 8.29. 

When Ze is below the critical value the front velocity 
relaxes quickly to a constant value, but when it is larger 
than Ze c the pulsations set in. This is illustrated in Figj4] 

This plot shows the result of integration of the follow- 
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0.0001 0.0002 0.0003 0.0004 0.0005 

x, cm 

FIG. 3: Coordinate dependence of different elements concen- 
trations for APROX13 nuclear network 



ing simplified system: 

dT_d_ 
dt dx 



dT 

dx 



RYl 



(13) 



Here, we assume the kinetic equation 
dt ~ RYa ' 



(14) 



and the rate of the reaction 7? = R(T) is taken to be 



R = cxp -E a 



1 



T Tn 



1 



(15) 



In our units the initial temperature is To, ashes have 
temperature TJ, = Tq + 1 and the activation energy is 



E a = Ze 



T b ~T 



= Ze(l+T ). 



(16) 



We have used the method of lines to simulate the equa- 
tion of thermal conduction, that is the space coordinate 
was discretized by finite differences, while integration 
over time was done by ODE integrators using methods 
of [49:] and [5(| ■ See details of the numerical techniques 
in more complicated situation of hydrodynamical evolu- 
tion in (5l| . Our simulations for this simple model give 
Ze c ~ 8.2. 

Timmes and Woosley Q simulated the flame fronts 
of different CO mixtures using one-dimensional time- 
dependent hydro code. They said nothing about this 
instability. At first glance this means that there is no 
TP instability in the degenerate presupernova condi- 
tions. However, later, Bychkov and Liberman [52[ and 
Nomoto et al. [53] have argued that C+C and O+O 
flame fronts should be TP-unstable. In the approxima- 
tion of delta function kinetics it is found in [52| that 
Ze c = 8 + 4-yo = 16.9, i.e. exactly by factor 2 larger 
than that for the Arrhenius law. It is even claimed by 



0.4 - 



0.3 



> 0.2 



0.1 




1000 2000 
time 



FIG. 4: Pulsating front velocity for Ze = 
obeying the Arrhenius law 



.3 and one reaction 



[52j that Timmes and Woosley [6J have overlooked the 
instability in their simulations. 

We believe that the latter extreme proposition should 
be rejected and we seek a real explanation of the apparent 
CO front stability taking into account that Ze for basic 
reactions of CO burning can indeed be higher than the 
critical value. The following study of stability is based in 
part on the analytical model by [54|. 

It is known that CO nuclear burning consists of the 
basic (C,a), (C,p), (0,a), ... reactions accompanied by 
secondary (p,7) and (a,j) reactions. It is very important 
that: 

1. Ze corresponding to the secondary reactions is 
rather low because of the lower height of the 
Coulomb barrier; 

2. the amount of energy yielded by the secondary re- 
actions is comparable to the yield of the basic ones. 

It is very difficult to investigate analytically any mod- 
els with many reactants. It should be made numerically 
as it was done in the Section IIIII and by Timmes and 
Woosley [6j. However, it is possible to develop a simple 
analytical model [54| with one deficient reactant which 
possesses in part the main properties of CO burning men- 
tioned above, and which can help us to understand the 
qualitative behaviour of the pulsating front. Moreover, 
it is not hard to develop a numerical model with only 
two reactions with different values of Ze numbers, one of 
them above and the other below the critical value. This 
model explains how the seeming contradiction between 
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the works [HI, [53[ , and the simulations in Ref . Q can be 
removed. 

We may compare the front structure of the simplified 
models with the front structures of C(50%)+O(50%) and 
O(60%)+Ne(30%)+Mg(10%) mixtures obtained in @, 
comparing, e.g., the ratio of the thickness of the con- 
ductive zone to the thickness of the reaction zone which 
reflects the effective value of the Ze number [47J. Then 
we find that the front structure, calculated in [6j is anal- 
ogous to the model front structure, when Ze belongs to 
the range 2-5. This range is noticeably below the crit- 
ical value Ze c . This fact suggests the explanation for 
the stability of CO flame front observed by Timmes and 
Woosley @. Moreover, it shows that the secondary reac- 
tions of CO burning are very important from the view- 
point of the TP stability and that the models considered 
in [52|, HH are oversimplified. 

Figs |5][7] illustrate, that for large enough Ze-number the 
instability can develop even when the system is stabilized 
by the secondary reactions. 

We use a simplified two-reactant system: 



dT 
~dt 



d_ 

dx 



dT 

dx 



-qiR<iY A Y B . (17) 



0.3 



0.2 



0.1 - 



2 reactions: Zel = 12.0 Ea2= 5.0 q2= 0.3 




200 400 
time 



FIG. 5: for Ze=12 



Here, we assume the following two kinetic equations: 
dY A 



dt 



dY B 
dt 



= -R 1 Yi-R 2 Y A Y B , 



R X Y%-R 2 Y A Y B , 



(18) 



(19) 



and the rates of the reactions Ri — Ri{T) are taken to 
be 



and 



R 2 = exp - — — 



E a 



T To + 1 



(20) 



(21) 



Examples in this section were given for the case of 
the switched-off hydrodynamics. Now we turn to the 
flames with an account of hydrodynamics, which is more 
relevant for realistic supernovae. 



V. FRONT PULSATIONS 

Here we will switch on hydrodynamics again. Accord- 
ing to the simulations in Section IIIII with the full set 
of reactions the instability behaviour was not observed 
(see the text below) . Now we will show that pulsational 
behaviour can exist in our simulations in the case of ar- 
tificially high Zeldovich numbers and will confirm that 
carbon and oxygen burning are stable. To carry out this 



task we should change the nuclear rate function. We will 
treat our reaction as one-step with rate as Arrhenius law: 



R{T) = Ae- B ' T *. 



(22) 



The constant B defines the Zeldovich number. Ze de- 
pends on the temperature of burning matter, according 
to definition (fj"2l): 



ZeAr 



B 



T 



9 burned 



(23) 



The value of constant A does not play any role here, but 
for definitencss we will determine it by the relation: 



T 2 T 2 
Ti Ti 



(24) 



A(B) = / R leal (T 9 )dT 9 



~ B / T °dT 9 



(25) 



By i?rcai(T) we assume here our variant with the simpli- 
fied nuclear network. The nuclear caloricity is preserved 
from this reaction too (we will demonstrate this explic- 
itly). From the results of the previous section we set the 
boundary temperatures as T\ = 1, T% = 10. Fig. [8] shows 
different rates used, compared to the "real rate" . We see 
that the rate with B w 50 fits the physical carbon burn- 



ing rate very nicely. This means that when T< 



9up 



10, 



8 




2000 




4000 6000 
time 



FIG. 6: for Ze=18 



FIG. 7: for Ze=24 



Ze w 5, which coincides with earlier estimates. Table 
IIIII presents results for different runs with po = 2 ■ 10 9 
g/cm 3 . We increase B from 20 to 200 and simultane- 
ously decrease caloricity q (from Ni to Mg) to avoid the 
flame acceleration. We see that for the normal caloricity 
(q = (5.6 — 9.2) • 10 17 erg/g) no pulsations appear for 
all variants of rates (upper part of the table). Zeldovich 
number, Ze, can be additionally increased by decreasing 
Tgbumed- This is done by caloricity change to q — 5 — 7 
MeV for one reaction (q = (2.0-2.8) ■ 10 17 erg/g). In this 
case pulsations appear for B > _B cr it — 112.5 (lower part 
of the table). Example of flame coordinate dependence 
is depicted in Fig. [9] 

So the result is that in "real" system (physical EOS 
and thermoconductivity) pulsations can exist, but only 
in the case of artificial parameters: the Zeldovich number 
should be increased by a factor of ~ 4. Table HVl presents 
the critical values of Ze number found in our numerical 
experiments for different densities of matter. 

In the case of complex nuclear network (many reac- 
tions) Ze number cannot be obtained in such straight- 
forward way, because it depends on the whole history of 
burning q(T). The front is stable in simulations with 
APROX13, presented in Section ITU for all runs. The only 
run must be considered carefully: That with initial 16 O 
and p = 2 • 10 s g/cm 3 demonstrates some irregular struc- 
ture in x(t) dependence (see Fig. [TTJl compare with Fig. 
[9]). This case is the best candidate for pulsations because 
of low Tf, and strong temperature dependence. But the 
period of these pulsations, n ~ 3 • 10~ 10 s, is close to 
T2 = Axpb/{u n p u ) = 8T0 -10 s, so it is probably a numer- 
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B=100.0 
B=125.0 
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FIG. 8: Comparison of artificial rates with 2 12 C— >Mg 2 



ical effect. Moreover, it does not contribute to physical 
effects like changing the flame velocity. 



VI. CONCLUSIONS 

The one-dimensional flame propagation in presuper- 
nova white dwarf has been considered. Flame proper- 
ties for different densities were obtained. It is shown 
that when only one nuclear reaction in the nuclear net- 
work is considered the flame velocity differs strongly from 
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2e-11 3e-11 4e-11 5e-11 6e-11 7e-11 8e-11 

t, s 

FIG. 9: Pulsational regime of flame with artificially increased 
Zeldovich number 



TABLE IV: Critical Zeldovich number for carbon burning at 
different densities: 



p, g 


/cm 3 


Ze cr 


2 


10 s 


18.4 < Ze < 21.4 


7 


10 8 


15.3 < Ze < 16.7 


2 


10 9 


13.7 < Ze < 14.1 



the more sophisticated net simulations. So this reaction 
should be used only in approximate simulations when the 
correct flame velocity is not required, or alternatively the 
nuclear rate should be revised to satisfy v n . 




0.0078 



6e-08 6.2e-08 6.4e-08 6.6e-08 6.8e-0i 
t 

FIG. 10: Flame coordinate for run with ie O, p = 210 8 g/cm 3 , 
APROX13 network 



We have studied also the one-dimensional pulsational 
instability. In "real" simulations no pulsation regime was 
obtained, so we used an artificial Arrhenius law to make 
it possible to change Zeldovich number, Ze. For high Ze 
numbers the pulsations do exist. By means of numerical 
simulations we have obtained critical Zeldovich numbers 
for densities in the range p — [2 -10 s , 2-10 9 ] g/cm 3 . These 
values are larger than those in real flame, which proves 
the one-dimensional stability of realistic flame fronts in 
supernovae. 



TABLE III: Results of simulations with Arrhenius law: 





A 


B q, 


10 17 erg/g Tgburned 


comm. 


1 


3.68 


10 6 


20.0 


9.2 


13.5 


flame 


2 


1.38 


10 8 


50.0 


9.2 


13.5 


flame 


3 


7.62 • 


10 12 


150.0 


5.6 


11 


flame 


4 


1.47- 


10 15 


200.0 


5.6 


14 


deton. 


5 


7.62 • 


10 12 


150.0 


2.8 


10 


puis. 


6 


3.60- 


10 10 


100.0 


2.8 


8 


flame 


7 


5.32 • 


10 11 


125.0 


2.8 


8 


puis. 


8 


1.39- 


10 11 


112.5 


2.8 


8.2 


flame 


9 


1.39- 


10 11 


112.5 


2.0 


7.8 


puis. 


10 


1.39- 


10 11 


112.5 


2.4 


8 


puis. 
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